Magnetic skyrmion lattices in heavy fermion superconductor UPt 3 
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Topological analysis of nearly S0(3) sp i„ symmetric Ginzburg-Landau theory, proposed for UPt3 by 
Machida et al, shows that there exists a new class of solutions carrying two units of magnetic flux: 
the magnetic skyrmion. These solutions do not have singular core like Abrikosov vortices and at 
low magnetic fields they become lighter for strongly type II superconductors. Magnetic skyrmions 
repel each other as 1/r at distances much larger then the magnetic penetration depth A, forming 
a relatively robust triangular lattice. The magnetic induction near H c \ is found to increase as 
(H - H cl ) 2 . This behavior agrees well with experiments. 

PACS numbers: 74.20.De, 74.25.Ha, 74.60. Ec, 74.70.Tx 



Heavy fermion superconductors have surprising prop- 
erties on both the microscopic and the macroscopic level. 
Charge carrier's pairing mechanism is unconventional. 
The vortex state is also rather different from that of s- 
wave superconductors: there exist unusual asymmetric 
vortices, and phase transitions between numerous vortex 
lattices take place [El . For the best studied material UPt3 
several phenomenological theories have been put forward 
[^|-Q which utilize a multicomponent order parameter. 
In particular, great effort has been made to qualitatively 
and quantitatively map the intricate H — T phase dia- 
gram. Most attention has been devoted to the region of 
magnetic fields near H c i- 

Magnetization curves of UPt3 near H c i are also rather 
unusual (see Fig.l). Theoretically, if the magnetization 
is due to penetration of vortices into a superconducting 
sample then one expects —AttM to drop with an infinite 
derivative at H c i (dotted line). On the other hand exper- 
imentally —AttM continues to increase smoothly (squares 
and triangles represent rescaled data taken from refs. ||] 
and |J] respectively). Such a behavior was attributed to 
strong flux pinning or surface effects Q. However both 
experimental curves in Fig.l, as well as the other ones 
found in literature, are close to each other if plotted in 
units of H c i . There might be a more fundamental expla- 
nation of the universal smooth magnetization curve near 
H c i. If one assumes that fluxons are of unconventional 
type for which interaction is long range then precisely 
this type of magnetization curve is obtained. 

Magnetization near H c i due to fluxons carrying TV 
units of flux <f>o = hc/2e, with line energy e and mu- 
tual interaction V(r), is found by minimizing the Gibbs 
energy of a very sparse triangular lattice: 
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where a = (<3?o/-Bv / 3) ^ is lattice spacing. When V(r) ~ 
exp[— Ar], the magnetic induction has the conventional 
behavior B ~ [log {H — H c x)] Q, while if it is long 
range, V(r) ~ then one finds B ~ (iJ — H c i) n+1 . 

The physical reason for this different behavior is very 



clear. For a short range repulsion, if one fluxon pene- 
trated the sample, many more can penetrate almost with 
no additional cost of energy. This leads to the infinite 
derivative of magnetization. On the other hand for a 
long range interaction making a place for each additional 
fluxon becomes energy consuming. Derivative of magne- 
tization thus becomes finite. 

It is generally assumed that although vortices in UPt3 
differ from the usual Abrikosov vortices in many details 
two important characteristics are preserved. First, 
their size A is well defined: magnetic field and interac- 
tions between vortices vanish exponentially beyond this 
length. Second, their energy is proportional to log/t. 
However, in this note we show on the basis of topological 
analysis of a model by Machida et al. || that there ex- 
ists an additional class of fluxons which we call magnetic 
skyrmions. They carry two units of magnetic flux N = 2 
and do not have a singular core, similar to ATC texture 
in superfluid 3 Hc |l(J. We show that their line energy 
e rs 2eq, £o = (<I>o/47rA) 2 is independent n and is smaller 
than that of Abrikosov vortices for strongly type II su- 
perconductors like UPt3 (k ~ 50). Magnetic skyrmion 
lattice becomes the ground state at low magnetic fields 
(H — H c i) / H c i *C 1. We further find that the repulsion of 
magnetic skyrmions is, in fact, long range: V(r) ~ 1/r. 
This allows us to produce a nice fit to the magnetization 
curve (solid line in Fig.l) which is universal (independent 
of k). 

The order parameter in the weak spin-orbit coupling 
model of UPt3 is a three dimensional complex vector: 
ijji{r) The Ginzburg-Landau free energy reads: 



AF, 



(2) 



F sym = -onkifi + ^(^iV*) 2 + f IV^I 2 (3) 
Ki (|2?^| 2 + \VM 2 ) + K 2 \V Z ^\ 2 + iB|, (4) 
AF = - 7 |^| 2 -A|V^| 2 + ^|^£M 2 (5) 



Ei= x ,y k{ (\V X ^ + \V y ^) + k\ \V Z ^ 



(6) 



where V 3 = d 3 - i(2e/hc)Aj and B 3 = (V X A)j. We 
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separated eq. (||) into a symmetric part F sym which is in- 
variant under the spin rotation group SO(3) sp j n acting 
on the index i of the order parameter, and into terms 
breaking the SO(3) sp i„ symmetry (anisotropy, coupling 
to antiferromagnetic spin fluctuations and spin-orbit cou- 
pling). Although AF is crucial in explaining the double 
superconducting phase transition in UPt3 at zero exter- 
nal magnetic field and the shape of H C 2(T) curve on the 
H — T phase diagram, it can be considered as a small per- 
turbation in the low temperature superconducting phase 
(phase B) well below its critical temperature T <C T~ ~ 
.45 K and at low magnetic fields H ~ H c i- Indeed, 
estimation of coefficients of AF at T = T c r/2 yields: 

7/a ~ .2, A/a ~ .05 and ^H 2 ^ / ^ 10~ 6 , and 

also k <C K Therefore, in a certain range of mag- 
netic fields and temperatures there is an approximate 
0(3) symmetry and we first turn to minimize F sym . 

In the vacuum of phase B the order parameter is 
if) = ipo(n + irn)/\f2,ip\ = a /Pi, ft _L rh, n 2 = rh 2 = 1. 
Stability requires: a > 0, (3i > and > 
The symmetry breaking pattern of phase B is as fol- 
lows. Both the spin rotations SO(3) sp i„ and the U(l) 
gauge symmetries are spontaneously broken, but a diag- 
onal subgroup U(l) survives. The subgroup consists of 
combined transformations: rotations by angle i9 around 
the axis I = H x fh accompanied by gauge transforma- 
tions e . Each vacuum state is specified by orientation 
of a triad of orthonormal vectors ft, rh, I. The vacuum 
manifold is therefore isomorphic to SO (3). Topological 
defects might be of two kinds: regular and "singular". 
To find regular topological line defects, it is enough to 
consider the London approximation pp| ], i.e. to assume 
that the order parameter gradually changes in space from 
one vacuum to another. The triad h, rh, I then be- 
comes a field. The Abrikosov vortex is a singular defect: 
it has a core where the modulus of the order param- 
eter vanishes and energy diverges logarithmically. Ac- 
cordingly, a cutoff parameter - the correlation length - 
should be introduced and one obtains log k dependence 
for the vortex line energy. This fact alone means that 
if there exists a regular solution it is bound to become 
energetically favorable for large enough k. Below we con- 
sider a situation when the external magnetic field is ori- 
ented along z axis and all configurations are translation- 
ally invariant in this direction. We use dimcnsionlcss 
units: r = Xr, A = ($ /2ttA) A, B = (<£>o/2ttA 2 ) B and 
F = (£o/2ttA 2 )F, where A = {^ q /2tt)-^J^koKI (the 
tilde will be dropped hereafter). The free energy takes 
form 

F L = 1/2 (dklf + (hd k m-A k ) 2 + B 2 (7) 
and the field equations are 

n p Vm p - A = Vx(vxl)=j, (8) 



At- f(f. Af) + 2j k (tx d k l) = 0. (9) 

Eq.(||) shows that the superconducting velocity is given 
by n p Vm p = — Vi?, where the angle r) specifies the ori- 
entation of vector h or rh in the plane perpendicular to 
I (see insert in Fig. 2). Thus, $ is the superconducting 
phase. 

Now we proceed to classify the boundary conditions. 
Magnetic field vanishes at infinity, while topology of the 
orientation of the triad n, rh, I at different distant points 
is described by the first homotopy group of vacuum man- 
ifold: 7Ti(SO(3))=Z2 Jl(|. It yields a classification of 
solutions into two topologically distinct classes ("odd" 
and "even"). This classification is too weak, however, 
for our purposes because it does not guarantee nontriv- 
ial flux penetrating the plane. We will see that con- 
figurations having both "parities" are of interest. In 
the presence of the magnetic flux possible configurations 
are further constrained by the flux quantization condi- 
tion. The vacuum manifold is naturally factored into 
SO(3)^SO(2)<g>S2 where the S2 is set of directions of 
I and the SO(2) is the superconducting phase For 
a given number of flux quanta N, the phase 1? makes N 
windings at infinity, see Fig. 2. The first homotopy group 
of this part is therefore fixed: 7i"i(SO(2))=Z. If, in addi- 
tion, I is constant, there is no way to avoid singularity 
in the phase 1? where \ip\ = 0. However, general require- 
ment that a solution has finite energy is much weaker. 
It tells us that the direction of / should be fixed only 
at infinity pl| . The relevant homotopy group is non- 
trivial: 7T2(S2)=Z. The second homotopy group appears 
because fixing / at infinity (say, up) effectively " compact- 
ifies" two dimensional physical space into S2. Unit vec- 
tor I winds towards the center of the texture. The new 
topological number is Q = (I/871") / £jj I (drf x djtj d 2 r. 
Therefore, all configurations fall into classes character- 
ized by the two integers N and Q. For regular so- 
lutions, however, these two numbers are not indepen- 
dent. Upon integrating the supercurrent equation, eq.(|^) 
along a remote contour and making use of the identity 
£pqslp(dil q ){djl s ) = {din p ){djm p )— (dim p )(djn p ), we ob- 
tain: Q = N/2. We call these regular solutions magnetic 
skyrmions. 

The lowest energy solution within the London approx- 
imation corresponds to N/2 = Q = — I (or A^/2 = Q = 
+1). We analyze a cylindrically symmetric situation and 
choose the triad n, rh, I in the form: 

I = e z cosO(p) + e p sinO(/9), 

n = (e z sin Q(p) — e p cos G(p)) sin^ + e v cos <p, (10) 
rh = [e z sin Q(p) — e p cos O(p)) cos ip — e v sin ip, 

where p and ip are polar coordinates and O = e z l. Bound- 
ary conditions are: O(0) = tt and 8(00) = 0. The vector 
potential is given by A = A{p)e ip . The general form of 
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such a configuration is shown in Fig. 2. The unit vec- 
tor I (solid arrows) flips its direction from up to down 
as it moves from infinity toward the origin. The phase 
$ (arrow inside small circles in Fig. 2) winds twice while 
completing an "infinitely remote" circle. If in cq. (0) 
only the first term were present we would deal with a 
standard SO(3) invariant nonlinear cr-model JO]. Being 
scale invariant, it possesses infinitely many pure skyrmion 
solutions <d s (p;5) = 2 arctan((5//?), which have the same 
energy equal to 2 (in units of £o) for any size S of a 
skyrmion. However, in the present case the structure of 
the order parameter is more complex and the above de- 
generacy is lifted by the second and third terms of eq. 

i- 

Below we make use of the functions O s (/0; 5) to explic- 
itly construct the variational configurations. We show 
that as size of these configurations increases the en- 
ergy is reduced to a value arbitrarily close to the ab- 
solute minimum of e ms — 2. Substituting eq. ([Toj) into 
eq.(Q) and integrating over the x — y plane we obtain 
the energy of the magnetic skyrmion in the form: e ms — 
£ s + £cur + e mag , where £ s = J pdp (9' 2 /2 + sin 2 0/2p 2 ) , 
£cur = J pdp [(1 + cos0)/p + A] 2 and e mag 
| pdp (A/ p + A') 2 . The first term e s is the same as in 
nonlinear cr-model without magnetic field. It is bound 
from below by 2, the energy of a pure skyrmion. The 
second term e cur , the " supercurrent" contribution, is pos- 
itive definite. One still can maintain zero value of this 
term when the field O(p) is a pure skyrmion s (p;<5) 
of certain size S. Assuming this one gets: A(p) = 
-(1 + cos6)//9 = -2p/(p 2 + 8 2 ). The third term, the 
magnetic field contribution (which is also positive def- 
inite), becomes e mag = 8/3<5 2 . It is clear that when 
5 — > oo we obtain energy arbitrarily close to the lower 
bound: e ms < 2 + 8/3<5 2 — * 2. Single magnetic skyrmion 
therefore blows up. 

If many magnetic skyrmions are present, then their 
interactions can stabilize the system. They repel each 
other, as we will see shortly, and therefore form a lattice. 
Since they are axially symmetric, the interaction is axi- 
ally symmetric and thus a triangular lattice is expected. 
Assume that the lattice spacing is a. At the boundaries of 
the hexagonal unit cells the angle O is zero, while at the 
centers it is n. The magnetic field B is continuous on the 
boundaries. Therefore, to analyze a magnetic skyrmion 
lattice we should solve eqs.(|8j)— (^) on the unit cell with 
these boundary conditions demanding that two units of 
flux pass through the cell (by adjusting the value of mag- 
netic field on the boundary). We have approximated the 
hexagonal unit cell by a circle of radius R = 3 3 / 4 a/v / 7r 
and the same area, and performed numerical integration 
of the equations 

... A' A l + cos6 , , 

A" + A =0, (11) 



n/l 6' sinG/2 + cose A „ 

6" + — + + 2A) = 0, 12 

p p V p J 

which follow from the cylindrically symmetric Ansatz of 
eq.©. Calculations for R from R = 5 till R = 600 
were done by means of a finite element method. The 
energy per unit cell in a wide range of R is satisfactory 
described (deviation at R = 10 is 1%) by the function 
e ce ii ~ 2 + 5.62/i?. Note that in the limit R — > oo we re- 
cover our previous variational estimate: e ce u — > e ms = 2. 
The dominant contribution to magnetic skyrmion energy 
at large R comes from the first term e s , similar to the 
analytical variational state described above. The contri- 
bution to e ce u from magnetic field, £ ma g, is small for large 
R but becomes significant in denser lattices. The most 
interesting feature of the solution is that the supercurrent 
contribution e cur to the energy of magnetic skyrmion is 
negligibly small for all considered values of R. This is 
to be compared with the usual Abrikosov vortex where 
at high k the total energy is dominated by magnetic and 
supercurrent contributions which arc of the same order 
of magnitude. Most of the flux goes through the region 
where the vector I is oriented upwards. In other words, 
the magnetic field is concentrated close to the center of 
a magnetic skyrmion. 

Line energy of Abrikosov vortices e v for the present 
model was calculated numerically (beyond London ap- 
proximation) in fl3| |. For k = 20 and 50 we obtain 
2e v /e ms 3.5 and 4.4 respectively. Therefore we expect 
that the lower critical field of UPt3 is determined by mag- 
netic skyrmions: h c \ = s ms /2N. Returning to physical 
units, 



H cl = $ /4^A 2 



(13) 



To find magnetization, we now utilize eq. ([!]). Interac- 
tions among magnetic skyrmions follows easily from the 
energy of a unit cell of the hexagonal lattice: V(r) = 
2(£ C eii — 2)/6 ~ 1.87/r. The resulting averaged magnetic 
induction, in units of $o/27rA 2 , reads 



B ~ 0.25 (H/H cl - 1) 



(14) 



This agrees very well with the experimental results, see 
Fig.l. For fields higher then several H c \ London approx- 
imation is not valid anymore since magnetic skyrmions 
will start to overlap. In this case one expects that ordi- 
nary Abrikosov vortices, which carry one unit of magnetic 
flux, become energetically favorable. The usual vortex 
picture has indeed been observed at high fields by Yaron 
et al. @ Curiously, our result is similar to conclusions of 
Burlachkov et al. |lq] who investigated stripe-like (quasi 
one dimensional) spin textures in triplet superconduc- 
tors. Having established the magnetic skyrmion solution 
of F sym we next estimated how it is influenced by var- 
ious terms of AF eqs.(^)-@. It was found that these 
perturbations do not lead to destabilization of a mag- 
netic skyrmion. 



3 



In conclusion, we have performed a topological classi- 
fication of the solutions in SO(3) sp j„ symmetric GL free 
energy. This model, with addition of very small symme- 
try breaking terms, describes heavy fermion supercon- 
ductor UPt3 and possibly other p-wave superconductors. 
A new class of topological solutions in weak magnetic 
field was identified. These solutions, magnetic skyrmions, 
do not have normal core. At small magnetic fields the 
magnetic skyrmions are lighter then Abrikosov vortices 
and therefore dominate the physics. Magnetic skyrmions 
repel each other as 1/r at distances much larger then 
magnetic penetration depth forming a relatively robust 
triangular lattice. H c \ is reduced by a factor log k as 
compared to that determined by usual Abrikosov vortex 
(see eq.((l|)). 

The following characteristic features, in addition to the 
slope of the magnetization curve, can allow experimental 
identification of a magnetic skyrmions lattice. 

1. Unit of flux quantization is 2$o- 

2. Superfluid density \ij)\ 2 is almost constant through- 
out the mixed state. This can be tested using STM tech- 
niques. 

3. Due to the fact that there is no normal core, in which 
usually dissipation and pinning take place, one expects 
that pinning effects are reduced. 

It is interesting to note that our results are actually 
applicable to another model of UPt3 with accidentally 
degenerate AE representations |ll| . Although this model 
adopts the strong spin-orbit coupling scheme, it has a 
structure closely related to F sym of eqs.(||)-(^) at low 
temperatures where both order parameters become of 
equal importance and can be viewed as a single three 
dimensional order parameter. 
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FIG. 1. Magnetization of magnetic skyrmion lattice (solid 
line) and experimental data for UPt3 from Ref.5 (squares) and 
Ref.6 (triangules) . Magnetization of Abrikosov vortex lattice 
(dotted line) is given for comparison. 




FIG. 2. Configuration of a magnetic skyrmion. Solid ar- 
rows represent I field while the "clocks" show that phase $ 
rotates twice as an infrtely remote contour is circumvented. 
Insert explains defintions of and G. 
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